subroutine more_nodes (BIGbis, nbis, mbis, n_regridbis, m_regridbis, NODE_CHANGEbis,Nonvalue, numsurf)
 
 implicit none 
 
 integer :: nbis, mbis, cv, tv, SPACE, n_regridbis, m_regridbis, i, j, numsurf, choice, NODE
 parameter (SPACE=6100)
 parameter (cv = 2)
 parameter (tv = 2)
 parameter (NODE =6100)
 real, allocatable :: BIGbis(:,:)
 real, allocatable :: NODE_CHANGEbis(:,:)
 real :: Nonvalue
 
 
 
 n_regridbis = (2*nbis)-1
 m_regridbis = (2*mbis)-1
  
 do i = 1, n_regridbis
   do j = 1, m_regridbis
   
   if (MOD(i,2)==1.and.MOD(j,2)==1) then; 
        if(BIGbis(((i/2)+1),(j/2)+1)-1<Nonvalue) then;
        NODE_CHANGEbis(i,j)=Nonvalue
        else;
      NODE_CHANGEbis(i,j) = BIGbis(((i/2)+1),((j/2)+1))
      end if 
   else;
      
        if (MOD((i+j),2)==0) then; 
            if(BIGbis(((i/2)+1),(j/2)+1)-1<Nonvalue.or.BIGbis((i/2)+1,(j/2))-1<Nonvalue.&
            &or.BIGbis((i/2)+1,(j/2))-1<Nonvalue.or.BIGbis((i/2),(j/2)+1)-1<Nonvalue) then;
            NODE_CHANGEbis(i,j)=Nonvalue
            else;
            NODE_CHANGEbis(i,j) = (BIGbis(((i/2)+1),(j/2)+1)+BIGbis((i/2)+1,(j/2))+BIGbis((i/2)+1,(j/2))+BIGbis((i/2),(j/2)+1))/4 
            end if       
        else;
           
           if (MOD(i,2)>0) then;
                if(BIGbis(((i/2)+1),(j/2)+1)-1<Nonvalue.or.BIGbis((i/2)+1,(j/2))-1<Nonvalue) then;
                NODE_CHANGEbis(i,j)=Nonvalue
                else; 
                NODE_CHANGEbis(i,j)=(BIGbis((i/2)+1,(j/2)+1)+BIGbis((i/2)+1,j/2))/2
                end if
           else;
                if(BIGbis((i/2)+1,(j/2)+1)-1<Nonvalue.or.BIGbis((i/2),(j/2)+1)-1<Nonvalue) then;
                NODE_CHANGEbis(i,j)=Nonvalue
                else;
                NODE_CHANGEbis(i,j)=(BIGbis((i/2)+1,(j/2)+1)+BIGbis((i/2),(j/2)+1))/2
                end if
           end if
        
        end if
        
    end if
    end do
  end do
  
  
 return
 
 end subroutine more_nodes